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1.  Introduction 


Predicting  the  moment  attributable  to  a  liquid  payload  in  a  spinning  and  coning  projectile  is  a 
problem  of  considerable  interest  to  the  Army.  Stewartson  (7)  was  first  to  consider  inviscid  pay- 
loads  contained  in  a  right  circular  cylinder,  and  his  results  gave  coning  Eigen  frequencies  that  can 
possibly  cause  catastrophic  yaw  for  a  liquid-carrying  projectile.  First  order  viscous  boundary  layer 
corrections  of  the  Stewartson  theory  were  made  by  Wedemeyer  (2)  and  Murphy  (3).  A  method  for 
calculating  the  linear  liquid  moment  using  the  full  linear  viscous  equations  with  boundary  layer 
corrections  confined  only  to  the  end  caps  (candlesticks)  was  presented  by  Hall,  Sedney,  and 
Gerber  ( 4 ,  5). 

A  further  interest  to  the  anny  is  to  consider  a  series  of  uniform  circular  cylinders  stacked  end  to  end 
separated  by  impenetrable  end  caps.  These  candlesticks  may  be  situated  along  the  symmetry  axis 
or  offset  from  this  axis  but  parallel  to  the  symmetry  axis  of  the  projectile.  Coning  motion-induced 
liquid  moments  are  considered  here  for  a  number  of  candlestick  configurations.  The  Eigen 
frequencies  for  such  configurations  are  shown  to  be  identical  to  those  found  by  Stewartson  (7). 

Liquid  payloads  contained  in  a  highly  permeable  material  have  been  of  interest  to  the  U.S.  Army 
for  some  time.  Laboratory  tests  and  flight  tests  have  shown  that  a  highly  permeable  medium  can 
significantly  reduce  the  spin-up  time  of  a  liquid  payload  (6,  7,  8).  Flight  stability  for  liquid- 
saturated  permeable  payloads  has  also  been  examined  by  D’Amico  ( 9 ,  10).  An  investigation  by 
Cooper  (77)  considered  inertial  waves  in  a  coning  projectile  with  a  saturated  porous  media  payload 
so  that  the  media  are  homogenous  and  isotropic.  The  present  work  extends  the  Stewartson  and 
Cooper  problems  by  considering  a  cylindrical  cavity  filled  with  a  permeable  medium  that  is 
impregnated  with  an  inviscid  liquid,  but  this  medium  is  not  necessarily  isotropic.  Following 
Cooper,  we  introduced  a  further  modification  by  segmenting  the  cavity  along  the  symmetry  axes 
into  a  sequence  of  equal  length  cylinders.  Each  of  these  cylinders  is  separated  by  impermeable  end 
caps.  The  porous  media  are  modeled  by  a  drag  term,  which  is  proportional  to  the  liquid  velocity 
relative  to  the  assumed  ridge  porous  media  that  are  added  to  the  linearized  Euler  equations.  This 
analysis  examines  the  induced  liquid  moment  as  a  function  of  parameters  found  by  Stewartson  (7) 
plus  parameters  describing  the  porous  media  and  the  number  of  segments  in  the  cylindrical  cavity. 


2.  Equations  of  Motion  for  the  Off-Axis  Candlestick  Configurations 


Figure  1  shows  the  X',  Y',Z'  axes  rotating  unifonnly  about  X'  with  angular  velocity  P  =  (P,  0,  0), 
and  figure  2  presents  details  of  the  internal  configuration  of  the  candlestick  payloads.  The  liquid  is 
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assumed  to  be  initially  rotating  as  a  rigid  body  with  the  same  angular  speed  P  so  the  velocity  V'  of 
the  liquid  inside  the  cylinders  is 

V'  =  Px(x,  R0cosB  +  rcos9,  R0 sin B  +  r sin 0)  (1) 

The  unperturbed  state  for  equation  1  satisfies  the  Euler  equation: 

p 

P2exx(exx(x,  RocosB  +  rcos0,  RosinB  +  rsin0))  = -V—  (2) 

P 

for  which  Ps  is  the  unperturbed  liquid  pressure.  Following  Stewartson  and  letting  the  position 
vector  have  the  Cartesian  form  R  =  (x,  R0  cos  B  +  r  cos  0,  R0  sin  B  +  r  sin  0)  allows  equation  2  to 
be  integrated  as  a  scalar  potential  so  that 

Ps/p=P2Rocos(B-0)-P2r2/2  (3) 
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Figure  2.  Details  of  payload  configuration. 


Now  we  excite  wave  motion  in  the  rigidly  rotating  liquid  by  perturbing  the  angular  velocity  of 
the  projectile  for  small  angular  components,  coy,coz,  so  that  the  total  angular  velocity  takes  the 

form 

Q  =  P  +  (0,  coy,  toz).  (4) 

The  projectile  is  assumed  to  undergo  small  angle  coning  motion  about  the  X',  Y',Z'  frame  and  is 
related  to  the  projectile  body  axes,  X,  Y,Z,  by  the  following  transfonnation  (3): 
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where  the  coning  damping  rate  is  s  the  coning  frequency  is  T  and  K0  is  the  magnitude  of  the 
small  coning  angle.  Equation  5  shows  that  the  angular  velocity,  to  first  order  in  K0 ,  is  written  as 
a  column  vector  with  body  frame  components 


P 

-  KO  P  T  [s  sin(t  P  (T  - 1))  +  cos(t  P  (T  - 1))]  eePTt 

-  KO  P  T  [sin(t  P  (T  - 1))  -  s  cos(t  P  (T  - 1))]  eePTt 


(6) 
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This  suggests  that  the  velocity  of  the  liquid  payload  should  be  given  by 

V  =  Px(x,  R0cosB  +  rcos9,  RosinB  +  rsin0)  +  v  (7) 

for  which  the  components  of  v  and  the  pressure  perturbation  p  all  have  the  magnitudes  of  order  K0 . 


The  liquid  payload  is  assumed  to  have  low  viscosity  and  the  magnitude  P  is  assumed  to  be  large 
so  that  the  fluid  motion  is  adequately  described  by  Euler’s  equations.  (The  M864  projectile  in 
free  flight  with  a  water  payload  has  a  Reynolds  number  on  the  order  of  106 .)  Neglecting  all 
higher  order  terms  in  small  K0  makes  the  Euler  equations  take  the  following  fonn: 


dV 

dt 


2PW  +  Nf 

p  dr 
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dW 
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+  2PV  + 


1  dp 
pr  d0 


dU  1  dp 

dt  p  dx 


0 
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(8) 


in  which  (V,  W,  Uj  are  cylindrical  components  of  the  perturbed  velocity  for  the  perturbed 
pressure,  p . 


Normal  boundary  conditions  at  the  solid  wall  satisfy 

(v,  W,  u)*n  =  flxRs»n  (9) 

in  which  n  is  an  outward  unit  vector  on  the  wall  and  Rs  a  point  on  any  cylinder  wall.  Equation  6 
when  substituted  into  equation  9  causes  the  nonnal  boundary  conditions  to  become 

V  =  tHe(x  K0PSePSt_iPt_ie) 

U  =  -9?e( K0 P (elB R0  +  r eI0)s el>Sl  ">l  '") .  (10) 

i  =  tJ-1,  S  =  (s  +  i)T 

These  show  that  separable  solutions  to  equation  8  are  obtained  by 

V 
w 
u 

_p. 

and  solving  for  the  velocity  components  yields 


v 

w 

u 

p 


,pst-iPT-ie 


(ID 


4 


dp 


dp 


v  =  - 


r  (S  -  i)  +  2 
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dp 
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U~”pP(S-i) 

Using  the  continuity  equation,  V  •  (v,  w,  u)  =  0 ,  produces  the  following  equation  for  the 
pressure  p : 

2  d2p  dp  d2p  2  2d2p  „ 
r  — -+r  —  +  — —  -  r  cr — -  =  0 


(12) 


dr2  dr  d  02 
(S  —  3  i)(S  +  i) 


dx~ 


(13) 


a  =  - 


(s-it 


At  this  point  in  the  analysis,  it  is  useful  to  consider  each  cylindrical  candlestick  ranging  from 
-  C  <  x  <  C  to  consist  of  an  end-to-end  sequence  of  N  equal  length,  A  ,  cylinders  with 
impenetrable  end  caps  so  that  A  =  2C/N .  Applying  equation  10  to  each  of  the  sub-cylinders 

gives  the  following  Fourier-Bessel  series  for  the  pressure  in  the  given  candlestick: 


P 

A 


=  K0  P2  a2  p 


e  10  ^Akcos(7rk(C  +  x)/A)j1(7rk 

keodd 

+  (rxE/a2  +rD/a)e10  +  Fx 
1 6  (—  l)n  C  S2  ( S  —  3  i) 


n2  a  k  N  (jj(zz)(S  +  i)  -  zz  J0(zz)(S  -  i)) 


zz  =  n  a  ka  N/2  C 
E  =  (S  -  i)S 

D  =  2  C  (N  -  2  n  -  1)S2/n  a 
F  =  eiB  R0(S-i)s/a 


0  <  n  <  N, 


(14) 


3.  Candlestick(s)  Liquid  Moments 


The  moment  induced  by  the  liquid  contained  in  the  segmented  cavity  is  calculated  from  the  time 
derivative  of  the  angular  momentum  field.  Non-dimens ionalizing  the  moment  with  2  n  p  a4  C  P2 
makes  it  convenient  to  write  the  side  moment  components,  MY  ,MZ  ,  as  (3): 
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(15) 


MY  J  +MZK  =  TClm  (2  n  p  a4  C  P2)ko  ePSMPt 

CLM  =  clsm(t,  n,  c)  j  +  clim(t,  n,  c)k 


Therefore,  the  induced  moment  on  the  entire  candlestick  consisting  of  N  end-to-end  sub-cylinders 
of  length  A  is  given  by 


TCLM(2;ipa4CP2)=-y> 

n=0 


JR; 
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dW 
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dQ 


(16) 


Q  =  27ta~  A 


Note  that  symmetry  attributable  to  the  geometry  of  each  candlestick  and  equation  14  causes  the 
axial  component  of  CLM  to  integrate  identically  to  zero. 


The  calculation  of  equation  16  is  tedious.  The  details  are  not  given  here  (3,  11),  but  the  result  is 
the  following  expression: 


T  (My  +  i  Mz) 


iS 


"(l2  cos  B  eiB  N2  R02  +  i  (4  C2  +  3  a2)N2  - 
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12a2  N2 
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n  a2  N2  (S  +  i)  keodd  [j,(zz)(S  +  i) -  zz  J0(zz)(S  -  i)]k4 


zz  =  7rakaN/2C, 


(17) 


The  summation  terms  of  equation  17  show  that  large  over-turning  moments  can  result  if 

Jj(zz)(S  +  i)  -  zz  J0(zz)(S  -  i)  =  0 
zz  =  7takaN/2C,k  e  odd 


(18) 


This  is  the  same  criterion  first  found  by  Stewartson  ( 1 )  for  the  symmetrically  located  liquid 
payload,  i.e.,  at  R  =  0  .  This  says  that  the  Stewartson  tables  can  be  used  to  determine  such  Eigen 
frequencies  where  the  aspect  ratio  is  taken  to  be  2  C/a  N  . 

One  further  comment  regarding  equation  17  is  that  CLM  approaches  the  value  produced  by  a 
frozen  liquid  as  N  becomes  large,  provided  that  Eigen  frequencies  are  sufficiently  removed 
from  the  region  of  interest.  In  fact,  the  rate  of  approach  to  the  frozen  limit  is  very  rapid  since  it 
goes  as  l/N2 .  Typical  coning  frequencies  of  the  M864  are  presented  in  figure  3.  Figures  3  and  4 
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give  examples  of  moment  coefficients,  CLSM  ,CLIM ,  for  the  M864  projectile.  The  value  for  the 
axes’  offset  is  taken  to  be  R  =  2a  and  the  orientation  angle  is  B  =  tt/3  .  In  these  plots,  the  sharp 
spikes  correspond  to  the  increased  over-turning  moment,  CLSM  ,  that  is  encountered  when  the 

projectile’s  coning  frequency  passes  a  Stewartson  Eigen  frequency.  These  results  clearly  show 
that  even  for  small  values  of  N ,  the  moment  coefficients  are  very  close  to  the  frozen  limit 
values.  The  reason  for  this  is  the  aspect  ratio  for  these  examples  is  relatively  large.  Some 
mathematical  manipulations  of  equation  17  show  that  for  large  aspect  ratios,  the  values 
Clsm  >  CLim  approach  their  frozen  limit  values  for  all  values  of  N  as  long  as  equation  1 8  is  not 

satisfied.  Figure  3  shows  that  coning  rates  are  in  the  range  of  0  <  T  <  0. 1,  and  figures  4  and  5 
indicate  that  large  over-turning  moments  caused  by  the  liquid  payload  should  not  result  in  flight 
instabilities  for  the  parameters  examined  here.  However,  if  the  parameters  and  particularly,  the 
aspect  ratio,  were  to  change,  instabilities  could  result  (see  figure  6  as  an  example). 
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Figure  3.  Free  flight  coning  frequencies  of  the  M864  projectile. 
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Aspect  Ratio  (C/a)=23  Ro/a=2 
M864 


T 


Figure  4.  Values  of  CLSM  showing  Eigen  frequencies  as  a  function  of  T  and  N  . 
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However,  if  the  parameters  and  particularly  the  aspect  ratio  were  to  change,  instabilities  could 
result.  Figure  6  is  one  such  example. 


Aspect  Ratio  (C/a)=23  R0/a=2 
M864 


Figure  6.  Values  of  CLSM  showing  a  possible  existence  of  an  instability  attributable  to 
an  Eigen  frequency  <0.1. 


4.  Equations  of  Motion  for  the  Symmetry  Axis  Porous  Media  Configuration 


For  this  problem,  the  moment  arm  R0  =  0  so  the  position  vector  becomes  R  =  (x,  r  cos0,  rsind) 
since  the  candlestick  axis  is  the  symmetry  axis  of  the  projectile.  When  the  liquid  payload  axis  is 
coincidental  with  the  projectile  symmetry  axis,  the  analysis  is  simplified  if  we  move  to  the 
X',  Y',  Z'  frame  rather  than  the  body  frame  used  before.  This  modifies  the  transformation  given 

by  equation  5  so  that 
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and  X,  Y,  Z  are  coordinates  in  the  non-rolling  reference  frame.  The  velocity  of  the  media  obtained 
from  the  rotation  kinematics  is  now  given  by 
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v;  =  -xK0P(S-i)ePSt_ie 
V'=Pr  +  ixKOP(S-i)e 
V  -rKyPfS-ije1 


PSt-iG 


(20) 


pst-ie 


and  r,  0,  x  are  cylindrical  coordinates  in  the  non-rolling  frame.  The  Euler  equations  are  now 
modified  to  account  for  an  inviscid  liquid  flowing  through  porous  media.  This  modification 
assumes  that  flow  in  porous  media  can  be  represented  by  additional  terms  that  are  proportional  to 
the  liquid  velocity  relative  to  the  media  which  is  taken  to  be  rigidly  attached  to  the  coning 
projectile.  Using  the  factor,  ePSl  ,  for  each  independent  variable  allows  the  modified  Euler 
equations  now  written  in  cylindrical  components  to  have  the  following  form: 


(S  +  Ct  -  i)v  -  2  w  +  x  Ct  P  (S  -  i)  +  — —  =  0 

pP  dr 


(S  +  Ct  -  i)w  +  2  v  -  i  x  Ct  P(S  -  i) ^  =  0 

pPr 

-rCxP(S-i)  +  — ^-  =  0 
pP  dx 


(21) 


(S  +  Cx  -  i)u 


In  this  equation,  v ,  w ,  u  are  the  r,  0,  x  components  of  the  fluid  velocity  and  p  is  the  perturbation 
pressure.  Note  the  terms  that  are  proportional  to  the  constants  Ct ,  Cx  represent  drag  that  the 

porous  media  exerts  on  flow  passing  through  these  media.  Motivation  for  this  model  stems  from 
Darcy’s  Law  (12)  which  says  that  the  drag  force  caused  by  porous  media,  Dr ,  is  given  by 

Dr  =  -  —  Vr  =  -pP2  aCr 

k  Pa 

p  — »  dynamic  viscosity 

k  — >  porosity  (dimensions  of  length  (2)) 

Vr  — »  velocity  of  fluid  relative  to  porous  media 

p  — >  liquid  density 

Dr  — >  pressure  gradient  induced  by  the  porous  media  on  the  fluid  flow 
Cr  — >  dimensionless  constant  which  is  a  measure  of  the  pressure  gradient 


so  that  Cr  = 


P 


PkP 


Solving  equation  2 1  leads  to  the  velocity  field  expressions: 
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V  =  - 


2ip  — —  (S  — i)r  , 

dr  V  x  Ct  P(S  -  i) 


w  = 


rpP(S-3i)(S  +  i)  S  +  Ct+i 

i(S  -  i)p  +  2  r  —  , 

dr  ix  Ct  P(S  -  i) 


u  =  - 


rpP(S-3i)(S  +  i)  S  +  Ct+i 
dp 

dx  ,  rCxP(S-i) 


-  + 


pP(S-i)  S  +  C,-i 


(22) 


Invoking  the  fact  that  the  flowing  liquid  is  incompressible  causes  the  continuity  equation  to  give 
the  following  expression  for  the  perturbation  pressure  p 


r  — -+r— -p-r  a  — —  = 
dr2  dr  dx2 

_2_  (S  +  Ct  -  3  i)(S  +  Ct  +  i) 

(S  +  Ct  -  i)  (S  +  Cx  -  i) 


(23) 


Equations  13  and  23  show  when  Ct  =  0,CX  =  0  causes  a2  =  a2 .  Using  the  demands  of  equation  20 
modifies  the  coefficients  for  the  Fourier-Bessel  series  as  follows: 


16(-lJICS(S-i)(S  +  Ct-3i) _ 

j  7T  a  k2  N  (j,(zz)(s  +  Ct  +  i)-  zz  J0(zz)(S  +  Ct  -  i)) 
zz  =  7rakaN/2C 

E  =  -(S-i)= 

D  =  2C(-N  +  2n  +  l)(S-i)S/Na 
F  =  0 


(24) 


0  <  n  <  N 


5.  Porous  Media  Liquid  Moments 


The  procedure  (equation  16)  for  calculating  liquid  side  moments  gives  the  following  expression 
for  liquid-saturated  porous  media  in  N  cylinders,  separated  by  impermeable  end  caps,  located 
along  the  projectile  symmetry  axis: 
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(25) 


T(Clsm  +  iCLIM) 


_  4  C2  (S  - i)((N2  -  l)s2  +  (Ct  +  2 i)N2  S  +  (i Ct  -  l)N2)  | 
3a2N2(S+Ct+i) 
-12C2(s2-l)-4C2  +  3a2  (S-i)2 
12  a2  4 


128  C2  S2(S-i)(S+Ct-3i) 
n4  N2  a2(S+Ct+i) 


I 

keodd 


_ Ji(zz) _ 

[jj(zz)(S  +Ct+i)  -  zz  J0(zz)(S  +Ct-i)]k4 


zz  =  7tkaN/2C 


Taking  the  limits  of  this  expression  for  large  N  or  large  Ct  gives  the  frozen  limit  values  of 
CLSM  and  CLIM  given  by  (3) 


C 


LSM  F 


4  C2/a2  +  3 


1-T 


C, 


2 

1 

=  —  + 


1  4C2/a2  +  3 
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T(s2-l) 


(26) 


whenever  T  is  not  in  the  neighborhood  of  an  Eigen  frequency  S  =  (i  +  s)T  which  by  equation  25 
satisfies 

J^zzXs  +Ct+i)-zz  J0(zz)(s  +C-i)  =  0  (27) 


The  definition  of  a  given  in  equation  23  shows  that  the  Stewartson  tables  (1,  3)  can  again  be 
used  to  find  Eigen  frequencies,  provided  that  Cx  =  Ct  and  S  is  replaced  by  S  +  Ct  while  these 

tables  are  used.  In  general,  these  frequencies  S  will  have  complex  values,  s  *  0  ,  thus  indicating 
that  damping/undamping  will  occur.  The  cases  when  Cx  ^  Ct  also  result  in  s  ^  0  but  numerical 

methods  need  to  be  used  to  find  such  Eigen  frequencies. 

Scheidegger  (14)  has  tabulated  permeability’s  k  =  1.0x1  CT5  cm2 ,  and  assuming  that  the  media 
studied  here  for  the  M864  have  similar  penneabilities,  which  average  Ct ,  Cx  =  0.6 .  Sample  plots 
showing  that  CLSM  and  CLIM  again  approach  the  values  for  a  frozen  liquid  as  Ct  becomes  large 
(see  figures  7  and  8),  which  are  in  agreement  with  equation  26  when  s  =  0 . 


Typical  plots  showing  CLSM  changing  with  increasing  values  of  N  are  given  in  figure  9.  For 
this  plot,  the  vertical  axis  assumes  values  defined  as  the  difference  between  CLSM  and  the  frozen 
values  given  by  equation  26: 


AC  =C  -C 

^  '“  LSM  —  '“LSM  '“LSM  F 


(28) 


The  sharp  peaks  indicate  Eigen  frequencies  which  are  now  complex  since  Cx  =  Ct  =  1/3  and 
therefore  s  *  0  .  Evidently,  CLSM  approaches  the  values  CLSM  F  as  N  gets  larger  for  values  of 
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T  sufficiently  removed  from  any  Eigen  frequencies.  In  general,  such  Eigen  frequencies  are 
found  with  a  numerical  search  routine  in  the  complex  plane. 


Aspect  Ratio  (C/a)=23  (R0/a)=0  N=2  8=0 
M864 


Figure  7.  Values  of  CLSM  as  a  function  of  T  for  increasing  Ct  and  Cx  . 


Aspect  Ratio  (C/a)=23  (R0/a)=0  N=2  8=0 

M864 


Figure  8.  Values  of  CLIM  function  of  T  for  increasing  Ct  and  Cx  . 


13 


Figure  9.  Values  of  ACLSM  showing  possible  flight  instabilities  (Eigen  frequencies)  as  a  function 

of  T  .  (As  N  increases,  the  payload  behaves  like  a  frozen  [liquid]  payload  when  T  is 
far  from  Eigen  frequencies.) 


6.  Calculation  Method 


The  equations  of  the  last  sections  need  to  be  calculated  for  a  wide  range  of  flight,  geometry,  and 
porous  media  parameters  all  of  which  require  values  of  Bessel  functions.  For  small  values  of  |zz| , 

simply  using  power  series  expansions  of  each  Bessel  function  works  very  well.  Bessel  functions 
at  large  values  of  |zz|  were  obtained  by  asymptotic  expansions  (13).  Generally,  calculating  Bessel 

functions  for  complex  arguments  for  intermediate  values  of  |zz|  is  a  non-trivial  problem  and  the 

methods  used  here  employ  Gaussian  continued  fractions.  This  author  has  judged  that  a  further 
discussion  of  these  methods  is  not  appropriate  for  this  article  but  the  reader  should  be  aware  of  the 
numerical  difficulties  associated  with  calculating  complex  Bessel  functions. 


7.  Conclusions 


The  off-axis  candlestick  problem  has  been  shown  to  be  equivalent  to  the  inviscid  Stewartson  (1) 
problem  whenever  porous  media  are  not  present  or  can  be  ignored.  Resonant  frequencies  are 
independent  of  the  candlestick  off-axis  position  and  can  be  found  with  tabulations  that  have 
already  been  found  ( 1 ,  3).  Values  not  found  in  such  tables  can  readily  be  obtained  by  simple 
numerical  root-finding  methods.  The  design  configurations  examined  here  show  that  flight 
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instabilities  should  not  occur  for  candlesticks  with  a  low  viscosity  liquid  and  as  N  increases,  the 
liquid  behaves  more  like  a  frozen  liquid. 

Cases  when  the  candlestick  contains  saturated  porous  media  that  are  located  along  the  symmetry 
axis  of  the  projectile  can  possibly  force  resonances.  In  these  cases,  the  Eigen  S  frequencies 
assume  complex  values  so  that  damping  or  un-damping  can  occur.  For  the  particular  case  when 
CT  =  Cx ,  it  is  possible  to  find  the  resonances  from  the  Stewartson  tables,  but  if  CT  ^  Cx ,  then  a 

numerical  search  for  resonances  in  the  complex  plane  for  the  Eigen  frequencies  is  required.  The 
current  analysis  is  based  on  assumed  values  of  the  porosities,  CT  Cx ,  and  the  results  indicate 

that  a  flight  instability  of  the  M864  is  unlikely  to  transpire.  However,  if  the  actual  porosity  of  the 
proposed  media  differs  significantly  from  the  assumptions  used  here,  catastrophic  flight 
instabilities  could  result.  A  thorough  search  for  “problem  Eigen  frequencies”  requires 
experimental  measurements  of  CT  and  Cx  in  order  to  gain  a  better  idea  where  to  search,  in  the 

complex  plane,  for  such  frequencies  (7,  8,  10).  In  all  cases,  the  liquid  moments  approach  the 
values  for  a  frozen  liquid  with  increasing  values  of  N  if  Eigen  frequencies  are  not  present. 
Similar  results  also  apply  for  increasing  values  of  CT  . 


15 


8.  References 


1.  Stewartson,  K.  on  the  Stability  of  a  Spinning  Top  Containing  Liquid.  Journal  of  Fluid 
Mechanics  1959,  5,  Part  4,  577-592. 

2.  Wedemeyer,  E.H.  Viscous  Correction  to  Stewartson ’s  Stability  Criterion',  BRL  Report 
1325;  US  Army  Ballistic  Research  Laboratory:  Aberdeen  Proving  Ground,  Maryland,  June 
1966.  (AD  489687) 

3.  Murphy,  C.H.  Angular  Motion  of  a  Spinning  Projectile  with  a  Viscous  Liquid  Payload', 
BRL-MR-3194;  U.S.  Anny  Ballistic  Research  Laboratory:  Aberdeen  Proving  Ground,  MD, 
August  1982.  (AD  A1 18676)  (See  also  Journal  of  Guidance,  Control,  and  Dynamics,  Vol.  6. 
pp.280-286,  July  -August  1983.) 

4.  Gerber,  N.;  Sedney,  R.  Moment  on  a  Liquid-Filled  Spinning  and  Nutating  Projectile:  Solid 
Body  Rotation',  BRL-TR-02470;  U.S.  Army  Ballistic  Research  Laboratory:  Aberdeen 
Proving  Ground,  MD,  February  1983.  (AD  A 125 3 32) 

5.  Hall,  P.;  Sedney,  R.;  Gerber,  N.  Fluid  Motion  in  a  Spinning,  Coning  Cylinder  via  Spatial 
Eigen  function  Expansion',  BRL  Technical  Report;  U.S.  Army  Ballistic  Research  Laboratory: 
Aberdeen  Proving  Ground,  MD,  August  1997  (AD  A 19075 8) 

6.  D’Amico,  W.P.;  Mark,  A.  The  Application  of  a  Highly  Permeable  Medium  to  Reduce  Spin- 
Up  Time  and  to  Stabilize  a  Liquid-Filled  Shell',  BRL-MR-0285 1 ;  U.S.  Army  Ballistic 
Research  Laboratory:  Aberdeen  Proving  Ground,  MD,  July  1978.  (AD  A058595) 

7.  Miller,  M.C.;  Molnar,  J.W.  Laboratory  Studies  to  Improve  Flight  Stability  of  the  M825 
Projectile',  CRDEC-TR-153;  U.S.  Army  Chemical  Research,  Development  and  Engineering 
Center:  Aberdeen  Proving  Ground,  MD,  June  1990. 

8.  Miller,  M.C.;  Molnar,  J.W.  Laboratory  Flight  Stability  Evaluation  ofM825Al  Production 
Payloads',  CRDEC-TR-240;  US  Army  Chemical  Research,  Development  and  Engineering 
Center:  Aberdeen  Proving  Ground,  MD,  October  1990. 

9.  D’Amico,  William  P.;  Soencksen,  Keith  P.  Aeroballistic  Testing  of  the  M825  Projectile: 
Phase  VIFLarger  Radius  Felt  Wedge  Payloads',  ARBRL-MR-3586;  US  Ballistic  Research 
Laboratory:  Aberdeen  Proving  Ground,  Maryland,  April  1987. 

10.  D’Amico,  W.P.;  Clay,  W.H.  Flight  Tests  for  Prototype  Felt  Wedge/White  Phosphorous 
Improved  Smoke  Concept',  BRL-MR-02824;  U.S.  Army  Ballistic  Research  Laboratory: 
Aberdeen  Proving  Ground,  MD,  April  1978.  (AD  A054643) 


16 


1 1 .  Cooper,  G.  R.  Moment  Exerted  on  Coning  Projectile  by  a  Spinning  Liquid  in  a  Cylindrical 
Cavity  Containing  a  Porous  Medium ;  BRL-MR-3677;  U.S.  Army  Ballistic  Research 
Laboratory:  Aberdeen  Proving  Ground,  MD,  June  1988.  (AD  A195291) 

12.  Darcy,  H.  Les  Fontaines  Publiques  de  la  Ville  de  Dijon,  1886. 

13.  McLachlan,  N.W.  Bessel  functions  for  Engineers',  Oxford  University  Press,  London,  1955. 

14.  Scheidegger,  A.  E.  The  Physics  of  Flow  Through  Porous  Media',  University  of  Toronto 
Press,  Canada,  1960. 


17 


NO.  OF 
COPIES 

1 

(PDF 

ONLY) 

1 


1 

1 

1 

2 

1 

10 


ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


DEFENSE  TECHNICAL 
INFORMATION  CTR 
DTIC  OCA 

8725  JOHN  J  KINGMAN  RD 
STE  0944 

FORT  BELVOIR  VA  22060-6218 

US  ARMY  RSRCH  DEV  &  ENGRG  CMD 

SYSTEMS  OF  SYSTEMS 

INTEGRATION 

AMSRD  SS  T 

6000  6TH  ST  STE  100 

FORT  BELVOIR  VA  22060-5608 

INST  FOR  ADVNCD  TCHNLGY 
THE  UNIV  OF  TEXAS  AT  AUSTIN 
3925  W  BRAKER  LN  STE  400 
AUSTIN  TX  78759-5316 

DIRECTOR 

US  ARMY  RESEARCH  LAB 
IMNE  ALC  IMS 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARL  Cl  OK  TL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARL  Cl  OK  T 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

SCHOOL  OF  AEROSPACE  ENGINEERING 
GEORGIA  INST  OF  TECHNOLOGY 
ATTN  DR  M  COSTELLO 
ATLANTA  GA  30332 

CDR  ARDEC 

ATTN  AMSRD  AAR  AEM  C  K  CHUNG 
R  LEE  A  READDY  M  CORZO 
AMSRD  AAR  AEM  A  W  KOENIG 
B  WONG  (5  CYS) 

PICATINNY  ARSENAL  NJ  07806-5000 

MR  MILES  MILLER 
504  HAVERHILL  ROAD 
JOPPA  MD  21085-4319 


1  ROY  KLINE 

27  FREDON  GREENDELL  RD 
NEWTON  NJ  07860-5213 

ABERDEEN  PROVING  GROUND 

1  DIRECTOR 

US  ARMY  RSCH  LABORATORY 
ATTN  AMSRD  ARL  Cl  OK  (TECH  LIB) 
BLDG  4600 

2  COMMANDER 
US  ARMY  ECBC 

ATTN  AMSRD  ECB  RT  D  WEBER 
BLDG  E3516 
APG  EA 

1  DIRECTOR 

US  ARMY  RSCH  LABORATORY 
ATTN  AMSRD  ARL  WM  J  SMITH 
BLDG  4600 

2  DIR  USARL 

ATTN  AMSRD  AR  WM  B  M  ZOLTOSKI 
BLDG  4600 

1  DIR  USARL 

ATTN  AMSRD  AR  WM  BA  D  LYON 
BLDG  4600 

1  DIR  USARL 

ATTN  AMSRD  ARL  WM  BD  B  FORCH 
BLDG  4600 

4  DIR  USARL 

ATTN  AMSRD  ARL  WM  BC  P  PLOSTINS 
G  COOPER  B  GUIDOS 
P  WEINACHT 
BLDG  390 

1  DIR  USARL 

ATTN  AMSRD  ARL  WM  BD  MNUSCA 
BLDG  390 


18 


